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Based on the interplay between photoionization and Raman effect in gas-filled photonic crystal 
fibers, we propose a new optical device to control frequency-conversion of ultrashort pulses. By 
tuning the input-pulse energy, the output spectrum can be either down-converted, up-converted, or 
even frequency-shift compensated. For low input energies, Raman effect is dominant and leads to a 
redshift that increases linearly during propagation. For larger pulse energies, photoionization starts 
to take over the frequency conversion process, and induces a strong blueshift. We have found also 
that the fiber-output pressure can provide an additional degree of freedom to control the spectrum 
shift. 


Hollow-core (HC) photonic crystal fibers (PCFs) have 
been developed in the last decade to explore nonlinear 
light-matter interaction^ with a flexibility far-beyond the 
all-solid optical fibers P, 0- HC-PCFs with Kagome- 
style lattice have offered unbeatable advantages to host 
strong nonlinear interactions between light and different 
gaseous media Q. These microstructures are character¬ 
ized by having wide transparency range and a pressure- 
tunable zero dispersion wavelength (ZDW) in the visible 
range . The abundance of gases with different proper¬ 
ties have led to the observation of various novel nonlinear 
phenomena [S- 11|. For instance, efficient deep-ultraviolet 
radiation has been observed in argon-filled HC-PCFs Q, 
while ultrabroadband supercontinuum generation span¬ 
ning from the vacuum-UV to the mid infrared regime is 
obtained using hydrogen-filled HC-PCFs [l^ . 

Raman self-frequency redshift, which continuously 
downshifts the central frequency of a single pulse dur¬ 
ing pr opagat ion, has been first observed in solid-core 
fibers [isIfiH, and later also in molecular gases 0. 
Photoionization-induced self-frequency blueshift due to 
plasma generation have been predicted and demonstrated 


in HC-PCFs filled by argon gas [8|, l9j, lll| . The core of this 


study is to investigate the non-trivial interplay between 
the Raman and photoionization nonlinear effects in HC- 
PCFs filled by Raman-active gases, and demonstrate its 
potential in designing novel photonic devices. 

We consider the propagation of an ultrashort pulse in 
a HC-PCF filled by a Raman-active gas. The evolu¬ 
tion of the pulse electric field can be accurately described 
via solving the unidirectional pulse propagation equation 
(UPPE) that does not require the slowly envelope ap¬ 
proximation [I30, 


BF' (jj^ 


is the combined gas and fiber dispersion of the fun¬ 
damental mode (HEii), rig is the gas refractive index, 
d = (2.4048 c)^ /2r^, r is the effective-core radius of the 
fiber, Pnl = -P {Pnl {z,t)} is the total spectral nonlin¬ 
ear polarization, jP is the Fourier transform, t is time in a 
reference frame moving with Vg, Pnl = Pk + Pi + Pr, Pk, 
Pi, and Pr are the Kerr, ionization-induced, and Raman 
nonlinear polarization, respectively, Pk = and 

is the third-order nonlinear susceptibility. We will 
describe below Pi and Pr in more detail. 

The time-dependent photoionization-induced polariza- 
tion satisfies the following evolution equation [19| 
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where e, me are the electron charge and mass, Ui is the 
ionization energy of the gas under consideration, rie is 
the generated free electron density governed by the rate 
equation 


dtrie = W (z, t) {riT - rie), (3) 


with riT the total density of the gas atoms, and W is 
the ionization rate. Using the Ammosov-Delone-Krainov 
model. 


= I ^|2n-l exp [-C2/ |P|] . 

\P\ 


( 4 ) 


where ci = ujpC, (4a;p-\/2TOeP//e)^" tOp = Uj/h is the 
transition frequency, h is the reduced Planck constant, 
^ = 2^”/ [nF (n) F (n -I- 1)], n = ^JUhJUi is the effective 
principal quantum number, F is the gamma function, and 
Uh — 13.6 eV is the ionization energy of the hydrogen 
atom. 

The dynamics of a specific Raman transition I in gases 
can be determined by solviM the Bloch equations for an 
effective two-level system [20|, 


where z is the longitudinal propagation distance, to is 
the angular frequency, E{z,uj) is the spectral electric 
field, (3i = l/vg is the first-order dispersion coefficient, 
Vg is the group velocity, c is the vacuum speed of light, 
eo is the vacuum permittivity, /3 (w) = oj/c [rig — 
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where and pij are the elements of the 2x2 polariz¬ 
ability and density matrices, respectively, w = P 22 — P 11 is 
the population inversion between the excited and ground 
states, P 22 + Pii — 1 is the trace of the density ma¬ 
trix (which states the conservation of probability, as¬ 
suming long dephasing times, as is the case for many 
Raman gases), p 2 i = P 12 , ai 2 = <^ 21 , Ti and T 2 are 
the population and polarization relaxation times, respec¬ 
tively, and the Raman polarization is given by = 
[aiipii + a 22 P 22 + 2ai2Re (P 12 )] utE. The Raman ‘co¬ 
herence’ field is defined as Re(pi 2 )- For weak Raman 
excitation, pn « 1 and P 22 ~ 0, i.e. the second term in 
Pr can be neglected, while the first term increases the 
linear refractive index of the medium by a fixed amount. 
The total polarization induced by all the excitable Ra¬ 
man modes is Pr = X]; Er ■ 

The interplay between the ionization and Raman ef¬ 
fects can be investigated during the pulse evolution inside 
the fiber by solving the strongly nonlinear coupled Eqs. 
(misi) . We are interested in exploring this interplay when 
pumping in two different anomalous dispersion regimes: 
(i) near the ZDW; (ii) away from the ZDW. 


Pumping near the ZDW^ the soliton fission process is 
expected in this regime due to non-negligible values of 
high-order dispersion coefficients. Due to the self-phase 
modulation effect, Raman-induced nonlinear redshift and 
dispersive wave generation, the spectrum is broadened 
dramatically and a supercontinuum generation is ob¬ 
tained 2l|. Panels of Fig. [3] show the evolution of an 
ultrashort Gaussian pulse with full-width-half-maximum 
(FWHM) 25 fs and centered at 805 nm in a hydrogen- 
filled HC-PCF with a ZDW at 487 nm. Hydrogen have 
two Raman transition modes - a rotational mode with 
period 57 fs, and a vibrational mode with 8 fs 0. 

The temporal and spectral evolution in the absence of 
the ionization effects are depicted in Fig. |3Ka,b). Emitted 
solitons from the soliton-fission process will impulsively 
excite the rotational Raman mode, since their durations 
will be shorter than 57 fs. Hence, they will experience 
redshift and deceleration in the spectral and time do¬ 
main, respectively, where the soliton with the shortest 
duration shows the strongest deceleration and redshift. 
Also, these solitons will excite a lagging sinusoidal tem¬ 
poral modulation of the medium refractive index that 
can be thought as a ‘temporal periodic crystal’ 22|. The 
vibrational mode of hydrogen will result in an additional 
Kerr nonlinearity, since it has an oscillation period less 
than the duration of the solitons. In Fig. |3l(c), the to¬ 
tal Raman coherence field is shown, and is a superposi¬ 
tion of two sine waves of different oscillation frequencies, 
corresponding to the rotational and vibrational Raman 
frequencies. 


The dispersive wave is periodically ‘twisted’ and is ra¬ 
diated during propagation because of the strong non- 
instantaneous Raman interaction with its emitting soli¬ 


ton. The long relaxation time of Raman in gases (for in¬ 
stance, the rotational mode in hydrogen has T 2 = 0.43 ns 
at gas pressure 7 bar and room temperature [23| 1 allows 
the dispersive wave to feel the induced temporal crystal. 
In analogy with condensed-matter physics, this wave will 
behave similarly to a wavefunction of an electron in a 
solid crystal in the presence of a uniform electric field 
[ 2 ^ . The applied constant force on the temporal crys¬ 
tal arises from the soliton redshift that is accompanied 
with a uniform acceleration in the time domain. Hence, 
the twist and the radiation of the dispersive wave can 
be interpreted as a result of combined Bloch oscillations 
[ 2 ^ and Zener tunneling [ 2 ^ . In the deep-UV, there is an 
emission of a very-weak dispersive wave (faint line on the 
right of the spectrum) due to interaction between posi¬ 
tive and negative frequency components of the electric 
field M- 

Switching on the ionization nonlinearity does not 
change the qualitative picture of the spectrum evolution, 
as shown in panel (d). However, the spectral broadening 
is reduced due to the ionization loss that suppresses the 
pulse intensity due to plasma formation [^, Q. Due to 
photoionization effects, the very-weak dispersive-wave in 
the deep UV disappears, and another wave in the mid- 
infrared (faint line on the left of the spectrum) is emitted 
instead, due to the photoionization-induced dispersion 

M- 

Working in the deep-anomalous regime, i.e. far from 
the ZDW, allows the observation of a soliton that can 
have a clean self-frequency shift, without emitting a dis¬ 
persive wave that lies far-way from the soliton central 
wavelength. Hence, we expect that in this regime the in¬ 
terplay between photoionization and Raman effects can 
lead to interesting features. The evolution of a sech- 
pulse with an input energy 2.7 p3 and centered at 1.55 
pLui in a 25 cm long H 2 -filled PCF and a ZDW at 403 
nm is depicted in Fig. [2] Panels (a,b) show the tem¬ 
poral and spectral evolution of the pulse, while panels 
(c,d) show the temporal evolution of the Raman coher¬ 
ence and the ionization fraction. The net self-frequency 
shift depends on whether photoionization or Raman non¬ 
linearity is dominant. Photoionizaton-induced blueshift 
occurs only when pulse intensity exceeds the ionization 
threshold [9j. However, Raman-induced redshift takes 
place along the whole propagation with a rate that has 
an approximately linear-dependence on the pulse inten¬ 
sity (22l |. Initially, the pulse is intense enough to ionize 
the gas and generates enough electrons that contribute to 
the pulse blueshift. As the pulse shifts to shorter wave¬ 
lengths, the group velocity dispersion decreases and the 
Kerr nonlinearity increases. This results in an adiabatic 
pulse compression llj that increases the amount of ion¬ 
ization and blueshift. As clear from panel (b), photoion¬ 
ization is dominant during the first propagation stage 
along the fiber. However due to the concurrent ioniza¬ 
tion losses, the blueshift process ceases after a certain 
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FIG. 1. (Color online), (a) Temporal and (b) spectral evolu¬ 
tion of a Gaussian pulse propagating with central wavelength 
0.805 /rm, energy 3 /rJ and FWHM 25 fs in a Ha-filled HC- 
PCF with a gas pressure 7 bar, and a hexagonal core-diameter 
27 fim in absence of ionization effects, (c) Temporal evolu¬ 
tion of Raman coherence wave induced by the rotational and 
vibrational excitations. The inset shows the total Raman co¬ 
herence at 2 = 3.5 cm. (d) The effect of ionization on the 
spectral evolution of the Gaussian pulse, using the ionziation 
energy of the hydrogen molecule Ui = 15.4 eV [2^. Contour 
plots in (a), (b), and (d) are given in a logarithmic scale and 
truncated at -50 dB. White arrows point to the weak (b) UV 
negative, and (c) IR dispersive waves. The linear, nonlinear, 
and Raman parameters of hydrogen can be found in [l^ . 

propagation distance zion, allowing the redshifting pro¬ 
cess to be dominant in the rest of the fiber. Panels (a,c) 
depict how the pulse and the Raman-induced coherence 
switch between acceleration and deceleration. The fiber 
length is a key-element in determining the net frequency 
shift at the fiber output. 

To explore the strong interaction between Raman 
and photoionization effects in this regime, the energy- 
dependency of the output spectral profile of an initial 
sech-pulse centered at 1.55 ^m is portrayed in Fig. [3Ka). 
For small input energies, ionization threshold that al¬ 
lows for plasma generation is not reached. Hence, Ra¬ 
man nonlinearity is dominant; its vibrational mode mod¬ 
ifies the Kerr nonlinearity, while its rotational mode in¬ 
troduces linear self-frequency redshift. Around 1.75 fiJ 
photoionization-induced self-frequency blueshift starts to 
take place due to plasma formation. As the input pulse 
energy increases, the induced blueshift starts to com¬ 
pensate then take over the rotational Raman-induced 
redshift. Initially, the blueshift increases linearly with 
the input energy, then there is a slight decrease beyond 
about 2.4 /iJ. This is because the photoionization pro¬ 
cess at this latter range of energies occurs very close to 
the fiber input, allows the Raman process to almost be 
dominant along the whole fiber, resulting in enhancing 
the reshift of the spectrum. Increasing the input energy 


FIG. 2. (Color online), (a) Temporal and (b) spectral evolu¬ 
tion of a sech-pulse with central wavelength 1.55 /rm, energy 
2.7 /rJ and FWHM 30 fs in a Ha-filled HC-PGF with a gas 
pressure 7 bar, and an effectivel core-diameter 18 pm. Tem¬ 
poral evolution of the (c) Raman coherence waves induced by 
rotational and vibrational excitations, and (d) ionization frac¬ 
tion. Contour plots in (a) and (b) are given in a logarithmic 
scale and truncated at -30 dB. 


more and more, the pulse breaks up and the clear soli- 
ton self-frequency shift is suppressed. The unambiguous 
interplay between photoionization and Raman processes 
for small input energies in this case suggests the proposal 
of a novel tunable device that can be used for either 
frequency-down or frequency-up conversion. The pulse 
central frequency can be scanned over 400 nm range from 
about 1.3 - 1.7 pm by changing its input energy over the 
range 0.5 - 2.8 pJ. 

The ratio between the fiber length and zion is a key- 
element in determining the frequency-up/down conver¬ 
sion range. Increasing the fiber length will certainly in¬ 
crease the frequency-down conversion, since the Raman 
nonlinearity will become dominant eventually. The latter 
is because the ionization loss halts the blueshift after the 
first few centimeters of the fiber and allows the redshift 
process to rule in the remaining distance. Contrarily, 
decreasing the fiber length, the frequency-up conversion 
range is enhanced over the frequency-down conversion 
range. So, there is a trade-off for choosing the length of 
the proposed device. However, the flexibility of tuning 
the output pressure at the gas-filled fiber end offers an 
additional degree of freedom to overcome this problem 
. At equilibrium, the pressure distribution across the 
fiber is P (z) = [Pq -b j- (P/ — Pq )] where Pq and 
Pl are the input and output pressures, and L is the fiber 
length. Panel (b) of Fig. [3] shows the dependence of 
the final spectrum of an initial sech-pulse on the output 
pressure with a hxed input pressure and pulse energy. By 
increasing the output pressure, the spectrum is linearly 
shifted towards the red side. For instance, by tuning the 






















4 


3 


2.5 


2 


1.5 


0.5 


Spectral evolution 



1 1.5 2 


0 


-10 


-15 


-20 


-25 


-30 


A (/fm) 



1.2 1.4 1.6 1.8 2 2.2 


A (^tm) 


FIG. 3. (Color online), (a) Energy and (b) output-pressure 
dependencies of the final spectral profiles of different sech- 
pulses with central wavelength 1.55 pm, and FWHM 30 fs 
in H 2 -filled HC-PCFs with length 10 cm, input gas pressure 
7 bar, and an effectivel core-diameter 18 pm. (a) The input 
pulse energy is scanned from 0.5 to 3.2 pj, and the output 
pressure is fixed at 7 bar. (b) The output pressure is scanned 
from 1 to 15 bar, and the input pulse energy is fixed at 1.3 pJ,. 
Contour plots are given in a logarithmic scale and truncated 
at -30 dB. 


output pressure from 1-15 bar, the spectrum is red- 
shifted by about 250 nm. This allows to shorten the 
fiber length close to Zion to enhance the photoionization- 
induced blueshift. Then, the reduced redshift at lower in¬ 
put energies can be recouped by raising the output-fiber 
pressure. This enables the proposed device to have a rel¬ 
atively wide tunable range for both frequency-up/down 
conversions. 

In conclusion, we have investigated the interplay be¬ 
tween the nonlinear photoionization and Raman effects 
in gas-filled HC-PCFs in two different dispersion regimes, 
pumping near or away from the ZDW. In the former 
regime, photoionization tends to degrade the broadening 
of the supercontinuum generation due to the ionization 
loss. However in the latter regime, photoionization in¬ 
duces a clear self-frequency blueshift that works against 
the Raman-induced redshift. This strong interplay sug¬ 
gests the proposal of a novel device that can be used for 
either frequency-up/down conversion via tuning the in¬ 
put pulse energy and the gas-pressure at the fiber output. 
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